Improved percolation thresholds for rods in three-dimensional boxes 
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We improve our previous results for the percolation thresholds of isotropically oriented rods in three 
dimensional boxes. We prove again the applicability of the excluded volume rule in the slender-rod 
limit (radius/length— > 0). Other limits for the rod sizes are discussed and important finite-size 
effects are revealed. 
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We have previously revealed errors of earlier simu- 
lations 1^ for the percolation of permeable rods in three 
dimensions. We have also presented corrected results 
that confirm the excluded volume theory. In this 
short preprint we give more accurate numerical values for 
the percolation thresholds, obtained by extended com- 
puter simulations. Given the importance of such results 
in designing composite materials (electrically conducting 
polymers, fiber-inforced plastics) it is useful to summa- 
rize and disscus the obtained percolation thresholds. 

The studied system is formed by permeable sticks with 
the form of capped cylinders (cylinders of length L and 
radius R capped with two hemispheres of radius R). It 
was conjectured Q that the percolation threshold, qp, 
is proportional to the inverse of the expected excluded 
volume. Vex- 
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(We denoted by Nc the number of sticks at percolation 
and by V the volume of the cube in which the percola- 
tion problem is considered.) For sticks with the form of 
capped cylinders the excluded volume is given by 
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Vex = ^R^ + SttLR^ + 'iL^R< sin{j) >, (2) 



where < sin('y) > is the average value of sin('-f) for two 
randomly positioned sticks, and 7 is the angle between 
them. For isotropic orientation of rods (< szn(7) >= 
7r/4) it was shown by a cluster expansion method |^] 
that the proportionality in (1) becomes equality in the 
— > slender-rod limit. This equality is known as 
the excluded- volume rule. 

We have studied the problem inside a cube with unit 
sizes. In order to preserve homogeneity near the cube's 
frontiers, the coordinates of the centers of the cylinders 
were generated uniformly in the interval [—{L/2 + R),l + 
{L/2+R)]. The orientation of the cylinders was isotropic. 
We detected the intersection of two capped cylinders by 
calculating the minimum distance between points on the 



two axes of the corresponding cylinders and by compar- 
ing this distance with 2R. Each time a new stick was 
generated, it was assigned to a cluster if it intersected a 
previous one, otherwise a new cluster was created; a stick 
may also unify two or more clusters that were previously 
separated. We considered percolation achived when a 
cluster spanned the cube from one face to the opposite 
one. The critical concentration is given by the number 
of sticks inside the cube at percolation, Nc. If a capped 
cylinder was only partially inside the cube, it contributed 
to Ne with a fractional value less than one, correspond- 
ing to the fraction of its volume inside the cube to its 
total volume. We generated many percolation realiza- 
tions (usually 5000) for each considered L and R pair, 
and the average Nc was calculated as the mean of the 
obtained percolation thresholds. 

The obtained results are summarized in Figs. 1-3. 
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FIG. 1. s = QpVex — 1 as a function of the R/L aspect 
ratios of the sticks. Data for three different stick lengths, L, 
are presented. 

On Fig. 1 we plot the quantity s = qpVex — 1 as a 
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function of R/L for various fixed L values. In the limit 
R/L ^ our simulations suggest the analytically pre- 
dicted s — relation This is nicely proven by our 
large-scale simulation data for L — 0.2 (for systems of up 
to Nc = S-IO'^ rods). The convergence for the applicabil- 
ity of the excluded volume rule is however rather poor. 
As an example in this sense we found in the — > 
limit {R/L < 0.05), for L = 0.15, s scaling as a function 
of R/L with an exponent of 0.58. From Fig. 1 is also 
clear that for smaller values of L and the same R/L ra- 
tios the value of s gets smaller. There are thus important 
finite-size effects, which are less evident in the R/L 
limit. 

We have checked that in the limit of L ^ the s = 
equality still does not hold. This is clear from our 
simulation data for R/L — 0.5, R/L = 0.25 and R/L — 
0.1. The data presented on Fig. 2 suggest that in the 
limit L ^ 0, s is linearly converging to 1.58, 1.46 and 1.14 
for R/L = 0.5, R/L = 0.25 and R/L = 0.1, respectively. 
This linear convergence may be useful for approximating 
the s value given R and L. 



• R/L=0.5; slope=5.00; intercept=1 .58 
A R/L=0.25; slope=1 .80; inlercept=1 .46 

* R/L=0.1; slope=0.36; intercept=1.14 
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FIG. 2. Finite size effects: s — f{L) for two different 
aspect ratios of the sticks. Continuous lines are the best linear 
fit. In the L — * limit we got s — 1.58, s — 1.46 and s = 14 
for R/L = 0.5, R/L = 0.25 and R/L = 0.1, respectively. 
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FIG. 3. s — f{R/L) for two different radii of the sticks. 

In conclusion, our results suggest that the excluded- 
volume rule is applicable in the slender- rod limit {R/L — > 
0) for stick lengths much smaller than the characteristic 
box-size. For all other cases important deviations from 
this rule are present. Decreasing the rod length improves 
the applicability of the excluded-volume rule, but the 
R/L — » condition is still necessary in the L — » limit. 
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Problems with the applicability of the excluded-volume 
rule appears also in the L ^ \ limit, where the length 
of the sticks approaches the charactheristic length of the 
box. This is nicely visible from Fig. 3 where we have 
plotted s as a function of R/L for two fixed R values 
{R = 0.01 and = 0.02, respectively). In the limit of very 
small R/L values one can observe again the increasing 
difference from the s = theoretical result, due to L 
values approaching the characteristic box size. 
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